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Abstract 

Since the pioneering works of Peleg and Scliaffer [5^, Althofer et al. [31, and Chandra et al. [llj . 
it is known that for every weighted undirected n- vertex m-edge graph G = (V, E) , and every integer 
fc > 1, there exists a {2k — l)-spanner with ©(ni+i/*^) edges and weight 0{k ■ n(i+<^)/'=) • uj{MST{G)), 
for any e > 0. (Here uj{MST{G)) stands for the weight of the minimum spanning tree of G.) Nearly 
hnear time algorithms for constructing (2fc— l)-spanners with nearly 0(n^+^/'^) edges were devised in 
[9l I34 [ l33 ] . However, these algorithms fail to guarantee any meaningful upper bound on the weight of 
the constructed spanners. 

To our knowledge, there are only two known algorithms for constructing sparse and light spanners 
for general graphs. One of them is the greedy algorithm of Althofer et al. [4 , analyzed by Chandra et al. 
[TT] . The drawback of the greedy algorithm is that it requires 0(to- + n-logn)) time. The other 
algorithm is due to Awerbuch et al. 6,, from 1991. It constructs 0(fc)-spanners with 0{k ■ ni+i/fc . A) 
edges, weig ht 0(fc2 • nVfc . A) • uj{MST{G)), within time 0(m • k ■ n^l^ ■ A), where A is the logarithm 
of the aspect ratio of the graph. 

The running time of both these algorithms is unsatisfactory. Moreover, the usually faster algorithm 
of pays for the speedup by significantly increasing both the stretch, the sparsity, and the weight of 
the resulting spanner. 

In this paper we devise an efficient algorithm for constructing sparse and light spanners. Specifically, 
our algorithm constructs ((2fc — 1) • (1 + e))-spanners with Oik ■ n-'^+^Z*^) edges and weight 0(fc • n^/*^) • 
uj{M ST{G)), where e > is an arbitrarily small constant. The running time of our algorithm is 
0{k ■ m + min{n • logri, m • a(n)}). Moreover, by slightly increasing the running time we can reduce 
the other parameters. These results address an open problem from the ESA'04 paper by Roditty and 
Zwick i34j. 
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1 Introduction 



1.1 Centralized Algorithms. Given an undirected weighted graph G = {V,E) and a parameter 
t > 1, a subgraph H = {V,E') of G {E' C E) is called a t-spanner if for every edge e = {u,v) G E, 
distffiujv) < t ■ distG{u,v). (Here distG{u,v) stands for the distance between u and v in the graph 
G.) Graph spanners were introduced in 1989 by Peleg and Schaffer [29] and Peleg and Ullman [30] . 
who showed that for every unweighted n- vertex graph G = {V,E) and an integer parameter A: > 1, 
there exists an 0(fc)-spanner with 0{n^^^^'') edges. Althofer et al. [4J improved and generalized these 
results. They analyzed the natural greedy algorithm for constructing graph spanners, and showed that 
for every n-vertex weighted graph G = {V,E) and an integer parameter k>l, this algorithm constructs 
a {2k — l)-spanner with 0(71^+^/*^) edges. (This is near-optimal [29].) They also showed that the weight 
of the resulting spanner is 0{n/k) ■ uj{MST{G)). (We will use the normalized notion of weight, called 
lightness^ which is the ratio between the weight of the spanner and uj{MST{G)).) In SoCG'92, Chandra 
et al. [H] improved the lightness bound, and showed that spanners obtained by the greedy algorithm 
have lightness 0{k-n^^'^''^/^), for an arbitrarily small e > 0. However, the running time of their algorithm 
is 0(m • {v}^^/^ + n ■ logn)), where m stands for \E\. Around the same time Awerbuch et al. [6] devised 
an algorithm that constructs 0(/c)-spanners with 0{k ■ n^+^l^ ■ A) edges and lightness 0{k^ ■ v}^^ ■ A), 
where A is the logarithm of the aspect ratio of the input graph. The running time of the algorithm of [6] 
is 0{m ■ k ■ n}/^ ■ A). 

In the two decades that passed since the results of [291 EOl IH El E] graph spanners turned out to be 
extremely useful. Among their applications is compact routing [301 [311 [36] , distance oracles and labels 
[281 l35l [33] , network synchronization [5] , and computing almost shortest paths [l3l [Ml [13 [13 [20] . Graph 
spanners also became a subject of intensive research for their own sake [131 113 113 13 [SZl [23 [13 133 [33 [8] . 

In particular, a lot of research attention was devoted to devising efficient algorithms for constructing 
sparse spanners for weighted graphs. Cohen [13j devised a randomized algorithm for constructing {{2k — 
1) • (1 + e))-spanners with 0{k ■ n^^^^'^) edges. Her algorithm requires expected 0{m ■ in}^^ ■ k) time. 
Improving upon [13], Baswana and Sen [9] devised a randomized algorithm that constructs {2k — 1)- 
spanners with expected 0{k-Ti}~^^/^) edges, within expected 0{k-m) time. Roditty et al. [33j derandomized 
this algorithm without any loss in parameters, or in running time. Roditty and Zwick [34J devised a 
deterministic algorithm for constructing {2k — l)-spanners with 0{n^^'^/^) edges in 0{k ■ n^+^Z*^) time. 
In the discussion section at the end of their paper Roditty and Zwick [31] write: 

"Another interesting property of the (original) greedy algorithm, shown by 111^ . is that the total weight 
of the edges in the {2k — 1)- spanner that it constructs is at most 0{n^^^^'^/'^ ■uj{MST{G)))]^ for any e > 0. 
Unfortunately, this property no longer holds for the modified greedy algorithm. Again, it is an interesting 
open problem to obtain an efficient spanner construction algorithm that does have this property. " 

In the current paper we devise such a construction. Specifically, our algorithm constructs {{2k — 1) • 
(1 + e))-spanners with 0{k ■ n^~^^^^) edges and lightness 0{k ■ n^^^), and does so in time 0{k ■ m + min{n • 
logn,m • a{n)}), where a{n) is an inverse Ackermann function. In other words, the running time of our 
algorithm is near-optimal and is drastically better than the running time 0{m ■ (n^"^^/'^ + n ■ logn)) of 
[llj and than that of [6] {0{m ■ k ■ n^^^ ■ A)). We pay for this speed up by a small increase (by a factor of 
(1 + e)) in the stretch and a small increase (by a factor of k) in the number of edges. (This comparison is 
with [llj. Our algorithm strictly outperforms the algorithm of [6j.) We also have another variant of our 
algorithm with a slightly higher running time {0{k ■ n^+^/'^)), and with 0{n^^^f^) edges and lightness 
0{k-n^/^). 

Note that the relationship between the stretch and lightness in both our results is essentially the 
same as in the state-of-the-art bound [H]. Specifically, in our result the slack factor (1 + e) appears in 
the stretch, while in [11] it appears in the exponent of the lightness. The number of edges in our slower 
construction (that runs in 0{k ■ n^^^^^) time) is the same as in [llj . The faster variant of our algorithm 
(that runs in near-optimal time of 0{k ■ m + min{n • logn, m ■ a(n)}), pays for the speedup by increasing 

^Actually, the weight in [11] is 0{k ■ n(i+")/'= • uj{MST{G))). 



1 



results 


stretch 


number of edges 


lightness 


running time 


Althofer et al. g] 


2fc - 1 


0{n'+"'') 


0{n/k) 


0(m ■ n^+i/'') 


Chandra et al. 1111 


2fc - 1 


0{n'+'"'') 


0{k-n'^'+'>"') 


0(m- (n^+i/* + n- log n)) 


Awerbuch et al. [5] 


0{k) 


0{k ■ n^+i/*^ • A) 


0{k' ■n"'" ■ K) 


0(m • k ■ n^'" ■ A) 


Our faster construction 


{2k -1) ■{! + €) 




Oik-n"-'^) 


0{k • m + min{n • log n, m ■ a(n)}) 


Our slower construction 


{2k -1) ■{! + €) 


0{n^+^/'') 


0{k-'n}''') 





Table 1: A concise comparison of previous and our constructions of light spanners. All the constructions mentioned in 
this table are deterministic. Our results are indicated by bold fonts. 

the number of edges by a factor of k. See Table [T] for a concise comparison of our and previous results 
on light spanners. (The lightness of all other spanner constructions [13 \ \9\ [33 l I34j is unbounded.) 

Chandra et al. jTT] also showed that the greedy algorithm gives rise to a construction of O(log^n)- 
spanners with 0{n) edges and constant lightness, and to a construction of 0(log n)-spanners with 0{n) 
edges and lightness O(logn). The running time of these constructions is 0{m ■ n ■ logn). Our algorithm 
also constructs spanners with the same (up to constant factors) parameters. The running time required 
by our algorithm to construct these spanners is 0(n^ • logn). 

1.2 Streaming Algorithms. In the streaming model of computation the input graph G = {V,E) 
arrives as a "stream", i.e., the algorithm reads edges one after another. The algorithm is required to 
process edges efficiently, and to store only a limited amount of information. In the context of computing 
spanners the natural memory limitation is the size of the spanner. Multi-pass streaming algorithms also 
allow several (ideally, just a few) passes over the input stream. 

The streaming model of computation was introduced by Alon et al. |2j and by Feigenbaum et al. 
|21j . The study of graph problems in the streaming model was introduced by Feigenbaum et al. [20]. In 
particular, Feigenbaum et al. [20] devised a randomized one-pass streaming algorithm for computing a 
(2/c + l)-spanner with expected 0{k-logn-n^^^/^) edges, using 0(A; • log n-n^/'^) processing time-per-edge. 
This result was improved in [T7] , who devised a randomized one-pass streaming algorithm that computes 
{2k — l)-spanners with expected 0{k ■ n^"'"^/'^) edges, using 0(1) processing time-per-edge. See also [7]. 
Elkin and Zhang |19] devised a multi-pass streaming algorithm for constructing sparse (1 + e, /3)-spanners. 
The number of passes in their algorithm is 0(/3). 

To our knowledge, there are currently no efficient streaming algorithms for computing light spanners. 
We show that our algorithm can be implemented efficiently in the streaming model augmented with the 
sorting primitives (henceforth, augmented streaming model). This model, introduced by Aggarwal et al. 
[1] in FOCS'94, allows to have "sort passes". As a result of a sort pass, in consequent passes one can 
assume that the input stream that the algorithm reads is sorted. (See [Ij for the justification of this 
model. The authors in [1] argue that "streaming computations with an added sorting primitive are a 
natural and efficiently implementable class of massive data set computations".) 

The algorithm of Chandra et al. [11] can be viewed as an algorithm in this model. After the initial 
sorting pass, it requires one pass over the input stream. As a result it constructs a {2k — l)-spanner with 
0(n^+^/'^) edges and lightness 0{k-n^^^''^^^), for an arbitrarily small e > 0. The processing time-per-edge 
of this algorithm is, however, 0{n^~^^^^), i.e., prohibitively large. 

We show that a variant of our algorithm computes ((2A;— l)-(l+e))-spanners with expected 0{k-n^^^^'') 
edges and expected lightness 0{k'^ ■ n^/^). It performs two passes over the input stream, that follow an 
initial sorting pass. In the first pass the worst-case (resp., amortized) processing time-per-edge of our 
algorithm is 0( io'giog„ ) (resp., 0{a{n))). The processing time-per-edge of our algorithm in its second 
pass over the input stream is 0{1). 

1.3 Our Techniques. Our algorithm is based on a transformation, which given a black-box construction 
of sparse (possibly heavy) spanners with a certain stretch t, efficiently produces sparse and light spanners 
with roughly the same stretch. We use this transformation in conjunction with a number of known 
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algorithms that produce sparse spanners, but do not provide any bound on their hghtness. 

Our transformation generahzes a metric transformation from [11]. Specificany, the metric transfor- 
mation of [Tl] converts constructions of sparse spanners for metrics into constructions of sparse and hght 
spanners (for the same metric). The generalized transformation that we devise appUes to weighted not 
necessarily complete graphs. (Observe that a metric can be viewed as a complete weighted graph.) 

There are a number of technical difficulties that we overcome in our way to the generalized transfor- 
mation. Next, we briefly discuss one of them. The construction of [llj hierarchically partitions the point 
set of the input metric into clusters. Then it selects a representative point from each cluster, and invokes 
its input black-box construction of sparse spanners on the metric induced by the representatives. One 
can try to mimic this approach in graphs by replacing each missing metric edge between representatives 
by the shortest path between them. This approach, however, is doomed to failure, as the overall number 
of edges taken into the spanner in this way might be too large. To overcome this difficulty we carefully 
select representative edges which are inserted into a certain auxiliary graph. Then the black-box input 
construction is applied to the auxiliary graph. As a result we obtain a spanner of the auxiliary graph, 
which we call auxiliary spanner. This auxiliary spanner Q = (p(,£) is a graph over a new vertex set U, 
i.e., U is not a subset of the original vertex set V. Next, we "project" the auxiliary spanner Q onto the 
original graph, i.e., we translate edges of £ into edges of the original edge set E. This needs to be done 
carefully, to avoid blowing up the stretch and lightness. Also, it is crucial that this translation procedure 
will be efficient. Interestingly, we do not project vertices of Vl onto vertices of V, but rather edges of £ 
onto edges of E. In particular, for a vertex u gU and two edges (u, x), (u, y) G <S, they may be translated 
into two vertex-disjoint edges {u' , x'), (u", y') G E. As a result a path in Q does not translate into a path 
in G, but rather into a collection of possibly vertex-disjoint edges. We show that these edges can be 
carefully glued into a path. This gluing, however, comes at a price of slightly increasing the stretch. 

1.4 Related Work. The large body of work on constructing graph spanners efficiently was already 
discussed in Section 1.1. The problem of constructing light spanners efficiently was also studied in the 
context of geometric spanners. See [15] and [23], and the references therein. 

1.5 Organization. In Section 2 we present and analyze our algorithm in the centralized model of com- 
putation. The algorithm is described in Section 2.1, and its analysis appears in Section 2.2. In Section 3 
we present a few variants of our basic algorithm (from Section 2.1). In particular, the streaming variant 
of our algorithm is presented in Section 3.5. 

1.6 Preliminaries. We will use the following results as a black-box (we write n = \V\,m = \E\). 

Theorem 1.1 f24\l [unweighted graphs] For any unweighted graph G = {y,E) and any integer k > 1, a 
{2k — l)-spanner with 0{n^~^^/^) edges can be built in 0{m) time. 

Theorem 1.2 |3 [weighted graphs I] For any weighted graph G = {V^E) and any integer k > 1, a 
{2k — l)-spanner with 0{k ■ n^"'"^/'^) edges can be built in 0{k ■ m) time. 

Remark: The algorithm of [9] is randomized, but was later derandomized in [33]. Henceforth, the 
algorithm provided by Theorem 11.21 is deterministic. 

Theorem 1.3 134^ [weighted graphs II] For any weighted graph G = {y,E) and any integer k > 1, a 
{2k — l)-spanner with 0{n^^^/'^) edges can be built in 0{k ■ n^^'^^^) time. 

Theorem 1.4 |j7| / [integer-weighted graphs] For any integer- weighted graph G = {V^E) and any integer 
k > 1, a {2k — l)-spanner with expected 0{k ■ n^^^^'') edges can be built in 0{SORT{m)) time, where 
SORT{m) is the time needed to sort m integers. 

Remark: The algorithm of [17J is randomized: while the guarantees 2k — 1 and 0{SORT{m)) on the 
stretch and running time, respectively, are deterministic, the guarantee 0{k ■ n^^^/^) on the size of the 
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spanner is in expectation. The fastest known randomized [26j (respectively, deterministic [25]) algorithm 
for sorting m integers requires expected 0(m- yloglogn) (resp., worst-case O(m-loglogn)) running time. 

We will henceforth refer to the algorithms of Theorems 11.11 11.21 11.31 and 11.41 as Algorithms U nwtdSp, 
WtdSp, WtdSp2, and IntWtdSp, respectively. 

2 The Basic Construction 
2.1 The Algorithm 

In this section we devise an algorithm LightSp that builds a light spanner efficiently. 

Let G = (V,E) be an arbitrary weighted graph, with n = \V\,m = \E\. Let A: > 1 be an integer 
parameter that determines the stretch bound of the spanner. 

We start with building an MST (or an 0(l)-approximate MST) T = {V,Et) for G. Let Mt = 
(y, distT) be the (shortest-path) metric induced by T. We then compute the Hamiltonian path C = 
{vi,V2, ■ ■ ■ ,Vn) of Mt drawn by taking the preorder traversal of T. Observe that for each 1 < i < 
71 — 1, distc{vi^ViJ^i) = distT{vi,Vi-^-i). More generally, for any pair u,v £ V of vertices, distc{u,v) > 
distT{u,v). Define L = a;(£); it is well known ([Ti], ch. 36) that L < 2-uj{T), and so L = 0{u{MST{G))). 

Let 1 < /9 < 2 be some parameter to be determined later, and define I = [log^n]. We partition 
the edges of E into i + 1 edge sets. The first edge set Eq contains all edges with weight in the range 
Wo = (0,^]. For each 1 < j < £, the jth edge set Ej contains all edges with weight in the range 

Define Vq = V, and let no = |Vo| = n. We use Algorithm WtdSp to build a {2k — l)-spanner 
H^ = {Vo,E'q) for Go = {V,Eo). 

Algorithm LightSp proceeds in i iterations j = 1,2, ... ,i. 

1. First, we divide the path C into nj = ^ = intervals of length /Xj = ^ each, where g^j^— j- < q < k 
is some parameter to be determined later. (Notice that rij > n, for all 1 < j < log^g + 1.) These 
intervals induce a partition of V in the obviou^ way; denote these intervals and the correspond- 
ing vertex sets by . . . , -^j"''^ and vj^\ ^'^^^ ■ ■ ■ ' respectively. While computing this 

(i) 

partition of V, we will store the index i of the vertex set Vj to which any vertex v £ V belongs 

in some variable indj(v), for each 1 < i < nj. For each interval 1 < i < rij, we define a new 

"dummy" vertex rj*"* (i.e., not present in G) which will serve as the representative of the interval 

The vertices rj*^ (respectively, intervals if^), 1 < ^ < nj, will also be referred to as the j- 
level representatives (resp., j-level intervals). For any vertex v £ V , its j-level representative is 
given by rj{v) = j-j*"'^^^''^^. We also define a dummy vertex set Vj, which contains all the j-level 

representatives. Observe that iVjl < min{n,nj} = min|n, 

2. We then compute an edge set Ej over Vj as follows. First, we remove from the jth edge set Ej all 
edges that have both their endpoints in the same j-level interval to obtain the edge set Ej . For each 
e = {u,v) £ Ej, let e = {rj{u),rj{v)) be a dummy edge of weight w(e) between the corresponding 
j-level representatives, and note that rj{u) / rj{y); we say that e = r(e) is the representative edge 
of e, and that e = s(e) is the source edge of e. Define Ej = {e = r(e) | e £ Ej}, and let Gj = iVj, Ej) 
be the corresponding multi-graph; note that Gj does not contain self-loops. Next, we transform Gj 
into a simple graph Gj = {Vj,Ej) by removing all edges but the one of minimum weight, for every 
pair of incident vertices in Gj. During this process, we store with every edge e £ Ej its source edge 
e = s(e) £ Ej. 

3. We proceed to building a {2k — l)-spanner H'j = {Vj,Ej) for Gj = {Vj,Ej) using Algorithm WtdSp. 
vertex that "lies" on the boundary of two consecutive intervals can be assigned to either one of them arbitrarily. 
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4. Next, we replace each edge e' G E'^ by its source edge s(e'). (Note that e' G C Ej is a dummy 
edge that connects some two distinct j-level intervals; denote them by I and Moreover, its 
weight (-o{e') is the smallest weight of an £'j-edge which connects I and I'. Hence the source edge 
s(e') is the minimum weight Ej-edge which connects the intervals / and /'.) Denote the resulting 
edge set by E* and define H* = {Vj,E*y, note that E* C Ej C Ej. We wiU refer to H* as the 
j-level spanner. 

Finally, we define E* = EtU E^U Uj=i -^j > ^■^d return the graph H* = (V, E*) as our spanner. 
2.2 The Analysis 

In this section we analyze the properties of the constructed graph H*. 

2.2.1 Stretch. We start with analyzing the stretch of the graph H* = {V, E*). 

Define str{k,q) = {2k - 1) • (1 + ^) + ^. Next, we show that the stretch of H* is at most str(k,q). 
In other words, we prove that distH*{u^v) < str{k,q) ■ distciu^v), for any edge e = {u^v) G E. 

If e G -Bq) then we have (iisiff* (ti, i;) < distH'^{u,v) < {2k — 1) ■ distG{u,v) < str{k,q) ■ distG{u,v). 

We henceforth assume that e e Ej, for some I < j < £. The next lemma shows that distH*{u,v) < 
str{k,q) ■ uj{e). (This lemma in conjunction with the triangle inequality imply that distH*{u,v) < 
str{k,q) ■ distG{u,v), for any edge e = {u,v) G E, which provides the required stretch bound on H*.) 

Lemma 2.1 There is a path of length at most str{k, q) ■ uj{e) in Hj U T = {V,EjU Et) between u and v. 

Proof: Since e is in Ej, its weight is in Wj = [^j, p ■ ^,]. In particular, we have u{e) > ^j. 

For any pair x,y £V of vertices, let IlT{x,y) denote the path in the MST T between them. 

Suppose first that u and v belong to the same j-lcvcl interval. Note that the length of all j-level 
intervals is Hj = so distT{u,v) < distc{u,v) < Hence, the weight of the path IIt{u,v) 

between u and in T is at most 

J < ^ < (2fc-l)-w(e) < str{k,q)-u{e). 

(The one before last inequality holds as q > 2^^-) Thus IIt{u,v) is a path of the required length in 
Hj U T between u and v. 

We henceforth assume that u and v belong to distinct j-lcvcl intervals. Let e = {rj{u),rj{v)) G Ej be 
the respective edge in Ej, where rj{u) and rj{v) are the j-level representatives of u and v, respectively. 
By the construction of Ej, we have w(e) < oj{e). Let W{rj{u),rj{v)) = {rj{u) = Uq,u'i, ■■■ ,u'f^ = rj{v)) 
be a path of length at most (2A; — 1) • distQ ,{rj{u), rj{v)) in H'- between rj{u) and rj{v). We have 

u{Y\'{rj{u),rj{v))) < {2k - 1) • dist^^{rj{u),rj{v)) < {2k - 1) • w(e) < (2A: - 1) • a;(e). (1) 

Notice that the edges of H'j (and of H' {rj{u),rj{v))) are not taken as is to the graph Hj. That is, 
each edge = {u'^,u^_^i) G Hj is replaced by its source edge s(e^) = {ui,Ui+i) G E*. To connect the 
two endpoints u^ and u^_^_i of the original edge G E^ in Hj U T, we take the path Il{u^,u^_^_i) = 
IlT{u'^,Ui) o {ui,Ui+i) o Ilj'{ui-^-i,u'j^^^) obtained from the concatenation of the path n7-(u',nj), the edge 
s(e') = {ui,Ui^i), and the path IlT{ui+i,u'j^^^). Observe that u'- is the j-lcvcl representative of Ui, i.e., 

u'j = rj{ui). In particular, they belong to the same j-level interval, and so uj{IlT{u[,Ui)) < ^. Similarly, 

uj{IlT{ui+i,u^_^_i)) < ^. Also, since s(e') £ E* C Ej, it holds that ^^(e^ > S,j. Recall that the weight of 
the source edge s(e9 = (uj,tii+i) of e[ is equal to that of = {u^,u'^_^_i), and so oj{s{e[)) = a;(e^) > ^j. 
Consequently, the weight io{U{u^,u'^_^i)) of the path Il{u'^, u'^_^i) satisfies 

w(n(^ii, Wi+i)) = uj{UT{u'i, Ui)) + u{ui, Ui+i) + u;{UT{ui+i,u'i^i)) 

< |+c.(eO + | < u:{e'^-{l + iy (2) 
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Let Il*{rj{u),rj{v)) be the path in H* U T between rj{u) and rj{v) obtained from Il'{rj{u),rj{v)) 
by replacing each edge {u^,u'^_^^l) in it by the path Il{u^,u'j^_^_l)) as described above. (See Figure [TJ (I).) 
Equations ^ and ([2]) yield 

ojiU*ir,iu),r,{v))) < a;(n'(r,(n), r,(z;))) • (^1 + ^) < (2A: - 1) • a;(e) • (^1 + ^) . 





Figure 1: (I) An illustration of the path Il*(rj{u),rj{v)). Edges = iui,Ui^i) of the path W{rj(u),rj{v)) in i/j are depicted 
by dashed curved lines. The MST-paths IlT{ui,Ui) £ T between vertices u'^ and iti that belong to the same interval are depicted 
by thick curved lines. All the other edges (i.e., edges in H*) are depicted by thin straight lines. (II) An illustration of 

the path Il*(u,v). The vertices u and rj{u) (respectively, v and rj{v)) belong to the same j-level interval, denoted here by /„ 
(resp., J„); the lengths of these intervals (denoted here by \Iu\ and \Iv\) are both equal to — . 



Finally, let U*{u,v) = UT{u,rj{u)) o Il*{rj{u),rj{v)) o I[T{rj{v),v) be the path in H* U T between 
u and V obtained from the concatenation of the paths IlTiu,rj{u)), H* {rj{u),rj{v)) and nj'(rj(f ), f ). 
(See Figure [TJ (II) for an illustration.) Since u and rj{u) belong to the same j-level interval, we have 
a;(nT(M, rj(u))) < ^. Similarly, u}{IlT{rj{v),v)) < Recall that uj{e) > ^j. We conclude that 

uj{U*iu,v)) = oj{UT{u,rjiu))) + oj{U*{rjiu),rj{v))) + oj{UTirjiv),v)) 

< k + ^2k-l)-u{e)-(l + -)+^ < u{e)-({2k-l)-(l + -)+-) 

q \ qJ Q \ \ qj qJ 

= str{k, q) ■ uj{e). 

Hence n*(n, v) is a path of length at most str{k, q) ■ w(e) between u and v in H* U T, as required. | 
2.2.2 Number of Edges. The next lemma bounds the number of edges in the spanner H* = (V, E*). 



Lemma 2.2 \E*\=o(k- n^+^A . 




Proof: Recall that \Vj\ < min{n,nj} = min|n, < n, for 1 < j < i- For an index 

logpQ + 1 < J < ^, we have -^^^ < 1. Hence n > -^f^ = nj, and so \Vj\ < nj. By construction, 
E* = EtUE'qUIJ^-^E*. Clearly \Et\ =n-l. ByTheoremO \E'^\ = 0{k ■ \Vo\^+^/'') = 0(A: • n^+V^). 
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Also, it holds that \E*\ = \E'j\ = 0{k ■ jl/jf +^/^) = 0{k ■ n^+Vfe)^ for each 1 < j < £. Next, consider an 
index log^ q + 1 < j < i; in this case l^jl < rij = -3^. Since q < k, we have q^^^ = 0(1). It follows that 

Since p > 1, we have 

logpg+l<j<i ^'^ ^ i=0 ^ i=Q ^< ' ^ 

It follows that 

I 

j=l l<i<logpg+l logpg+l<i<€ 

= 0(k . n^+y^) + 0(log,g • k ■ ni+iA') + O [k ■ q ■ n^+i/^) • ^ (^^Y ' 

= o(^k-n'+'/>'-(^l + -^^^ 

(Note that for g > 1, it holds that log^ q = 0{-^).) | 

2.2.3 Weight. The next lemma estimates the weight of H* = {V,E*). 

Lemma 2.3 uj{E*) = O ■ n^l^ ■ ^) • uj{MST{G)). 

Proof: First, observe that the edge weights in Eq are bounded above by ^. Since \Eq\ = 0{k ■ n^^^^^), 
it follows that 

Lo(E') < 0(k ■ n^+^/h ■ - = Oik-n^l^-L). (4) 

n 

The edge weights in Ej are bounded above hy ■ p = ■ for each 1 < j < i. 
Consider first indices j in the range 1 < j < logp^ + 1. Since \Ej \ = 0{k ■ n^^^^'^), we have (jj{Ej) = 
0{k ■ n^+^l^) . {pi . i^) = 0{k ■ -n}!^ ■ pi ■ L). It follows that 



J2 ^(Ej) = 0(k-n^l^-Ly Y P' = 0(k-n^/^ -L-qy (5) 

i<i<logp?+i i<i<logpg+i 

Next, consider indices j with log^ q + l < j < L By \E*\ = O q - ■ {^ ^^l^/^ y ^ . Hence 



Observe that 



It follows that 
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By Equations (H), (P and ([6]), 



i<i<logpg+i logpq+i<j<e 
= uj{T) + 0{k ■ n^l^ ■L)^o(k- n^l^ ■ L ■ + O (k^ ■ q- n^'^ ■ L ■ -^—^ 

= o(^k^ -n^/^ --^ycviMSTiG)). 

(Since q > ^f}^-^ , we have k"^ ■ q = ft{k).) | 

2.2.4 Running Time. The next lemma bounds the running time of our construction. 

Lemma 2.4 The graph H* can be built within time 0{k ■ m + min{n • logn, m ■ a(n)}). 

Proof: The tree T can be built in time 0(min{m + n ■ log n,m • a(m, n)}) [12] . Computing the 
Hamiltonian path C of Mt, as well as its weight L = oj{C), take another 0(n) time. Having computed 
L, we can partition the edges of E to the i+1 edge sets Eq, Ei, . . . , Ei in 0{m) time in the obvious way. 
By Theorem 11.21 building the spanner Hq for Go = {V,Eq) requires 0{k ■ \Eo\) = 0{k ■ m) time. 
Next, we bound the running time of a single iteration j in the main loop, for 1 < j < i- 

1. For each 1 < j < ^, we allocate an array Aj of size rij. For each x G Vj, the entry Aj[x] will contain 
the linked list of x's neighbors in the multi-graph Gj (with the weights of the representative edges). 
For each edge e = {u,v) £ E we determine the index j £ [£] such that a;(e) G Wj (and so e E Ej). By 
the locations of u and v on the path C, we determine the j-level representatives rj{u) and rj{v) of u 
and V, respectively. We also determine whether the edge crosses between different j-level intervals. 
If it is the case (i.e., e G Ej), then we insert the edge into the array Aj. (In other words, we update 
the linked lists of both rj{u) and rj{v) in Aj.) In this way we form the multi-graphs Gi, G2, . . . , 
within 0(|-E|) = 0{m) time. To prune these multi-graphs into simple graphs Gi, G2, . . . , G^, we 
need to remove all edges but the one of minimum weight, for every pair of incident vertices in Gj. 
This can be done deterministically within 0{m) time and space. (Alternatively, one can apply one 
of the algorithms for constructing {2k — l)-spanners with 0{k ■ n^"^^/*^) edges for weighted graphs 
[U [Ml [IZ] directly on the multi- graphs Gi,G2, ■ ■ ■ ,Gi, without pruning them first. It is easy to 
verify that this does not affect their guarantees for stretch, size and running time.) 

2. By Theorem [Ol the time needed to build the {2k - l)-spanner H'j = {Vj,E'j) for Gj = {V,Ej) is 
0{k-\Ej\) = 0{k-\Ej\). Summing over all ^ iterations, the overall time is ^^^^^ 0( A;- |ii^j|) = 0{k-m). 

3. Replacing each edge of Ej by its source edge takes 0(1) time, thus the graph H* = {V,E*) is built 
in 0(|£^^|) = 0(|i?j|) time. Summing over all £ iterations, the overall time is X]j=i Od-Bjl) = 0{m). 

It follows that the total running time of the construction is 

0(min{m -|- n • log n, m • a(m, n)}) -|- 0(A; • m) = 0(A; • m -|- min{n • log n, m • a(n)}). | 

We remark that there are randomized algorithms for constructing MST within 0{m + n) time |27j . 
Employing one of them reduces the running time of our algorithm to just 0{k ■ m). 

2.2.5 Summary. We summarize the properties of the spanner H* = (y, E*) in the following theorem. 
(We substituted p = 2 to optimize the parameters of the construction.) 

Theorem 2.5 Let G = {V-iE) be a weighted graph, with n = \V\,m = \E\. For any integer k > 1 
and any number < q < k, a (^{2k — 1) • (1 + |) -|- |^ -spanner with O {k ■ n^+^Z'' ■ (1 + g)) edges and 
lightness O (fc^ • n^^^ ■ qj, can be built in 0{k ■ m + min{n • logn,m • a{n)}) time. 
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By substituting q = in Theorem 12.51 for some smah constant e > 0, we obtain: 

Corollary 2.6 Let G = {V,E) he a weighted graph, with n = \V\,m = \E\. For any integer k>l and 
any constant e > 0, a {{2k — 1) • (1 + e))-spanner with O {k ■ n^^"^/^^ edges and lightness O {k"^ ■ n^/^), 
can he huilt in 0{k ■ m + min{n • log n, m ■ a{n)}) time. 

Note also that by substituting q = Q{k/e), we get a {2k — 1 + e)-spanner with O {k"^ ■ n^^'^/^') edges 
and lightness O {k^ ■ n^^^^, within the same time. 

3 Variants of the Construction 

In this section we devise some variants of the basic construction of Section [2] by applying a small (but 
significant) modification to Algorithm LightSp. 

Notice that Algorithm LightSp employs Algorithm WtdSp as a black box for building (i) the spanner 
i^o = {yo,K) for the graph Go = {V,Eo), and (ii) the spanner H'j = {Vj,E'j) for the graph Gj = {Vj,Ej) 
in step 3 of the main loop, for each 1 < j < i- Instead of employing Algorithm WtdSp, we can employ 
some of the other black-box spanners that are summarized in Section 1.6. 

3.1 First Variant (Slightly Increasing the Stretch) 

There is a significant difference between the spanner Hq, and the spanners Hj for j > 1. While the aspect 
ratio (defined as the ratio between the maximum and minimum edge weights) of the graph Go may be 
unbounded, the aspect ratio of all the graphs Gj,l < j < i, is bounded above by p. This suggests that 
we may view the graphs Gj as unweighted, and will thus be able to use Algorithm UnwtdSp for building 
the spanners ILj for them. As a result, the stretch of each spanner H'j will increase by a factor of p. 

Denote by ?t^* the variant of LI* obtained by employing Algorithm UnwtdSp for building the spanners 
Hj, I < j < L (We still use Algorithm WtdSp to build the spanner H'q.) 

We will next analyze the properties of the resulting construction %* . This analysis is very similar to 
the analysis of the basic construction that is given in Section 2.2, hence we aim for conciseness. 

Stretch. We first show how to adapt the stretch analysis of Section 2.2.1 to the new construction %* . 
The proof of Lemma [2.11 carries through except for Equation ([T]) that needs to be changed. Observe that 
all weights of edges in Ej belong to Wj = [Cj,p ■ Cj]- Hence they differ from each other by at most a 
multiplicative factor of p. Note that IL'j is an unweighted {2k — l)-spanner for Gj, i.e., a subgraph of Gj 
that contains, for every edge e G Ej, a path with at most {2k — 1) edges between its endpoints. It is easy 
to see that Lij provides a {p ■ {2k — l))-spanner for Gj. Hence, instead of Equation ([1]), we will now have 

oj{U'{rj{u),rj{v))) < p ■ {2k - I) ■ dist^^{rj{u),rj{v)) < p ■ {2k - 1) ■ uj{e) < p ■ {2k - 1) ■ uj{e). 

As a result, the upper bound on the weight of the path U* {rj{u),rj{v)) will also increase by a factor of 
p, and we will have uj{Il* {rj{u),rj{v))) < p ■ {2k — 1) • uj{e) ■ ^1 + 1^ Consequently, we will get that 

u{U*{u, v)) < uj{e) ■ • {2k - 1) • (^1 + l) + l) . Hence the stretch of Ti* is p ■ {2k - 1) • (1 + |) + |. 

Number of Edges. Next, we show how to adapt the size analysis for the new construction LL* . 
The only change from the proof of Lemma [2.21 is that now the upper bound on \Ej\ = \Ej\, for 1 < j < i, 
is smaller by a factor of k. The reason is that we now use Algorithm UnwtdSp rather than Algorithm 
WtdSp to build the spanner Ll'j = {Vj,Ej), for 1 < j < £; compare Theorem 11.11 with Theorem 11.21 
Note that the bound on \E'q\ remains unchanged, though, as we still use Algorithm WtdSp to build the 
spanner LI'q. We will get 

e 

\E*\ = \Et\ + \E'o\ + J2\Ej\ = n - 1 + 0{k ■ n'+'/'') + Yl 1^11+ Yl '^H 

i=i i<i<iogpg+i iogpq+i<j<e 
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0{k . n^+V^) + o(log, q . n^+^l^) + O (q ■ n^+V^') . ^ 



i-1 



logp g+i<j<e 



O n^+^/^ ■ k + 



q 



p-i 



Weight. We now show how to adapt the weight analysis for the new construction 71*. 

The only change from the proof of Lemma 12.31 is that now the upper bound on u}{E*), 1 < j < i, is 

smaller by a factor of k. The bound on oo{E'q) remains unchanged (see Equation Hence, 

u{E*) = u{T) + u{E',)+ ^(^1)+ E ^(^1) 

l<j<logpg+l logpg+l<j<^ 

= u{T) + 0{k ■ n^/^ ■L)+0 (n^l^ ■L-(^+o{k-q- v}/^ ■ L ■ -^—^ 

= o(^-n^'^ ■(^ + -^^^-^^-uj{MST{G)). 

Running time. Clearly, the running time of %* is not higher than that of the basic construction. 
We summarize the properties of the resulting construction %* in the following theorem. 

Theorem 3.1 Let G = (y, E) he a weighted graph, with n = \ V\,m = \E\ . For any integer k>l, and any 
numbers 2k-i < q < k and 1 < p < 2, a (^p ■ {2k — 1) • (1 + |) + |^ -spanner with O ^n^'^^^'^ ■ (^k + 

edges and lightness O Ik ■ n)-!^ ' ( ^ + 1^ ) ) ? '-^'^ built in 0{k ■ m + min{n • logn, m ■ a{n)}) time. 



P 

By substituting q = 6(7), /> = 1 + ©(e) in Theorem 13. H for an arbitrary constant e > 0, we obtain: 

Corollary 3.2 LetG= {V,E) be a weighted graph, withn= \ V\,m= \E\. For any integer k > 1 and any 
small constant e > 0, a {{2k — 1) • (1 + e))-spanner with O {k ■ n^"*"^/*^) edges and lightness O {k ■ n^^^^, 
can be built in 0{k ■ m + min{n • log n, m ■ a{n)}) time. 

Corollary 13.21 in the particular case k = O(logn) gives rise to an 0(log n)-spanner with 0{n- logn) edges 
and lightness O(logn). The running time of this construction is 0{m ■ logn). 

3.2 Second Variant (Increasing the Running Time) 

Denote by T-L the variant of H* obtained by employing Algorithm U nwtdSp for building the spanners Hj , 
1 < J < ^1 and employing Algorithm WtdSp2 (due to Roditty and Zwick [34]) to build the spanner H'q. It 
is easy to see that the stretch bound of the resulting construction 'H is equal to that of 7i*. Also, the size 
and weight bounds of T-l are better than those of H* and Ti*, but the running time (which is dominated 
by the running time of Algorithm WtdSp2) is higher. The analysis of this variant is very similar to the 
analysis of the basic construction (in Section 2.2) and the analysis of the first variant (in Section [3.ip . 
and is thus omitted. 

We summarize the properties of the resulting construction T-L in the following theorem. 

Theorem 3.3 Let G = {V,E) be a weighted graph, with n = \V\. For any integer k > 1, and any 
numbers 2k-i < Q < k and 1 < p < 2, a (^p ■ {2k — 1) • (1 + |) + -spanner with O (^n^^^^^ ■ ^1 + ^5~r^^ 

edges and lightness O (^k • n^^^ ■ ^^ri^; can be built in 0{k • n^^^^^) time. 

By substituting q = Q{\),p = 1 + ©(e) in Theorem 13. 3| for an arbitrary constant e > 0, we obtain: 

Corollary 3.4 Let G = {y,E) be a weighted graph, with n = \V\. For any integer k > 1 and any small 
constant e > 0, a {{2k — 1) • (1 + e))-spanner with O (n^+^Z'^') edges and lightness O {k ■ n^^^^ , can be built 
in 0{k ■ n^+-'^/'^) time. 
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Note that the tradeoff between the stretch and lightness exhibited in Corollary 13.41 is the same as in 
(There the stretch is 2k — 1, and the lightness is 0{k ■ n^^^''^^^).) The number of edges is 0{n^~^^/'') in 
both results, but our running time is 0{k ■ n^^^^*^), instead of 0{m ■ (n^^^/'^' + n ■ logn)) in 
By substituting q = Q(l/k),p = 2 in Theorem 13.31 we obtain: 

Corollary 3.5 Let G = {V, E) be a weighted graph, with n = \V\. For any integer k > 1, an 0{k'^)- 
spanner with O (n^"''^/'^) edges and lightness O (n^/'^), can be built in 0{k ■ ri^^^l^) time. 

Corollary 13.51 in the particular case k = O(logn) gives rise to an 0(log^ n)-spanner with 0{n) edges and 
lightness 0(1). We can extend this result to get a general tradeoff between the three parameters by 
substituting k = 0(logn),p = 2 in Theorem 13.31 and taking 1 < ^ = | = O(logn). 

Corollary 3.6 Let G = {V, E) be a weighted graph, with n = \ V\. For any number 1 < I = Oilogn), an 
0(logn • i)-spanner with 0{n) edges and lightness O(^^y^), can be built in 0{n^ • logn) time. 

The tradeoff of Corollarv l3.6l was also given by Chandra et al. [11], but their running time is O(m-n-logn). 

3.3 Third Variant (Increasing the Running Time Some More) 

Denote by Ti the variant of H* obtained by employing Algorithm WtdSp2 for building all the spanners 
H'j, 1 < j < £, as well as the spanner Hq. It is easy to see that the stretch bound of the resulting 

construction H is equal to that of the basic construction H*. Also, the size and weight bounds of H are 
exactly the same as those of the second variant T-L, but the running time is slightly higher (by a factor of 
(1 + g^)) than that of Ti. The analysis of this variant is very similar to the above, and is thus omitted. 

We summarize the properties of the resulting construction Ti in the following theorem. (We substituted 
/9 = 2 to optimize the parameters of the construction.) 

Theorem 3.7 Let G = {V,E) be a weighted graph, with n = \V\. For any integer k > 1, and any 
number < q < k, a (^{2k — 1) • (1 + |) + -spanner with O [^n^^^^'' ■ (1 + g)) edges and lightness 
O [k ■ n}/^ ■ q), can be built in 0{k ■ n^+i/'^ . (1 + q^)) time. 

3.4 Fourth Variant (Integer- Weighted Graphs) 

The fourth variant of our construction applies to integer-weighted graphs only. Denote by T-Lmt the 
variant of H* obtained by employing Algorithm LntWtdSp for building all the black-box spanners, i.e., 
the spanners Hj, 1 < j < £, and the spanner Hq. The new construction achieves exactly the same bounds 
as the basic construction H* , except for the running time. The running time of this variant consists of 
two parts. Specifically, it is the time required to compute the MST and the time required to compute 
a spanner. On an integer-weighted graph the first task can be done in 0(m + n) time [22], while the 
second task requires 0{SORT{m)) time (see Theorem 1.4 in |17J). Note, however, that the construction 
becomes randomized. 

We summarize the properties of the new construction T-Lint in the following theorem. 

Theorem 3.8 LetG = (V^E) be an integer-weighted graph, withn = \ V\,m = \E\. For any integer k > 1 

and any number 2k-i < Q < k, a (^{2k ~ 1) " (1 + |) + |^ -spanner with expected O [k ■ rji+i/^ . (1 -)- qj^ 

edges and expected lightness O [k'^ ■ n^/^ • q), can be built in 0{SORT(m)) time. 

3.5 Fifth Variant (Spanners in the Streaming Model) 

In this section we analyze our algorithm in the augmented streaming model. Specifically, this is the model 
introduced by Aggarwal et al. [Ij, which allows sorting passes over the input. 

Our algorithm relies on a streaming algorithm for constructing sparse (but possibly heavy) spanners 
from [ITj, summarized in Theorem 13.91 below. We remark that Theorem 13.91 also applies to multi-graphs. 
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Theorem 3.9 |j7| / For any unweighted n-vertex graph G = {V^E) and any integer k'>l, there exists a 
one-pass streaming algorithm that computes a {2k — l)-spanner with 0{k ■ n^~^^/^) edges (expected). The 
processing time-per-edge of the algorithm is 0{1), and its space requirement is 0{k ■ n^^^/^) (expected). 
Also, in the augmented streaming model (if the algorithm accepts a sorted stream of edges as input) the 
algorithm produces {2k — l)-spanners with 0{k-n^^'^/^) edges (expected) for weighted graphs as well. The 
processing time-per-edge and space requirements of this algorithm are the same as in the unweighted case. 

Our algorithm wih run £ + 1 copies of the algorithm for weighted graphs from Theorem 13.91 We will 
denote these copies Aj, for < j < i. 

Our algorithm starts with a sorting pass over the stream of edges. After this pass, in the consecutive 
pass the algorithm reads edges in a non-decreasing order of weights. The objective of the first pass (after 
the sorting pass) is to compute an MST of the input graph. To accomplish this, the algorithm maintains 
a Union-Find data structure (see Ch. 21 in [H]). It is known [IHIE] that all operations can be performed 
in worst-case O( io'giogra ) ti™^, and in total 0{n + m ■ a{n)) time, using 0{n) space. 

As a result of the first pass, the MST T is computed. The Hamiltonian path C of Mt is computed 
between the passes. In addition, the algorithm maintains the location on C of every vertex v G V. It also 
initializes i -\- 1 arrays Aj of size nj each, < j < £. 

Then the algorithm performs the second pass over the sorted stream of edges. For each edge e = (n, v) 
that the algorithm reads in the second pass, the algorithm determines the index j such that u}{e) G Wj. 
Then it tests if the edge crosses between different j-level intervals, i.e., belongs to Ej. If it does not, this 
edge is skipped. Otherwise the algorithm passes the edge e to the j'th copy Aj of the weighted streaming 
algorithm from Theorem 13.91 The streaming algorithm Aj will ultimately produce the spanner for 

Gj = {Vj,Ej). (Here we follow the notation of Section [2Tl ) If the streaming algorithm Aj decides to 
insert e into the spanner Hj, it will also insert its source edge s{e) into the ultimate spanner H* = {V, E*). 
As a result the spanner H* will contain the union of the j-level spanners, for < j < £. The edge set of 
the MST will also be inserted into H* (either between the passes, or after the second pass). 

This completes the description of the algorithm. By Theorem 13.91 the (expected) space requirement 
in the second pass is 

Ej=o 0{k ■ n]^^^'') = 0{k ■ ni+i/^). The processing time-per-edge is 0(1). 
We summarize our streaming algorithm in the following theorem. 

Theorem 3.10 In the augmented streaming model, for any weighted graph G = {V, E), any integer k > 1, 
and any constant e > 0, our algorithm requires two passes after the sorting pass. It computes a {{2k — 1) • 
(1 -|- e))-spanner with expected 0{k ■ n^'^^/^) edges and expected lightness 0{k^ ■ n^/^). The expected space 
requirement is 0{k ■ n^^^'^^^). The worst-case processing time-per-edge of the first (respectively, second) 
pass is 0( io'g^g„ ) (resp., 0{1)). Moreover, the overall processing time of the first pass is 0{m ■ a{n)). 

As was mentioned in the introduction, the algorithm of [11] can be viewed as an algorithm in this model. 
It constructs {2k — l)-spanners with ©(n^"*"^/*^) edges and lightness 0{k ■ n^^"*""^^/*^). It requires one pass 
after the initial sorting pass, but its processing time-peg-edge is very large (specifically, it is 0{n^^^/^)). 
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